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Abstract 

We introduce creation and annihilation operators of pseudo-Hermitian 
fermions for two-level systems described by pseudo-Hermitian Hamilto- 
nian with real eigenvalues. This allows the generalization of the fermionic 
coherent states approach to such systems. Pseudo-fermionic coherent 
states are constructed as eigenstates of two pseudo-fermion annihilation 
operators. These coherent states form a bi-normal and bi-overcomplete 
system, and their evolution governed by the pseudo-Hermitian Hamilto- 
nian is temporally stable. In terms of the introduced pseudo-fermion op- 
erators the two-level system' Hamiltonian takes a factorized form similar 
to that of a harmonic oscillator. 



1 Introduction 

In the last years a great deal of interest has been devoted to the study of non- 
Hermitian Hamiltonians with real spectrum (see [1]- [17] and references therein). 
Bender and Boettcher were the first to touch on this issue [7], by introducing 
the notion of PT-symmetry for one-dimensional non-Hcrmitian Hamiltonian 
H v = p 2 + x 2 {ixY , {u > 0), that possesses real, positive and discrete spectrum. 
In [T51H5] the Bessis - Zinn Justin conjecture about the reality of the spectrum of 
the PT-symmetric Hamiltonian —d 2 /dx 2 — (ix) 2M for M > 1 has been proven. 
A criterion for the reality of the spectrum of non-Hermitian PT-symmetric 
Hamiltonians is provided in [17] . 



By definition a PT-symmetric Hamiltonians H satisfies the relation 

[if, PT ] = , (1) 

where P and T are the operators of parity and time-reversal transformations 
respectively. These are defined according to 

PxP = -x, PpP = TpT = -p, TUT = -il , (2) 

where x,p,l are the position, momentum, and identity operators respectively, 
acting on the Hilbert space, and i := y—l. 

Later, Mostafazadeh [Ej-[H| introduced the notion of pseudo-Hermiticity in 
order to establish the mathematical relation with the notion of PT symmetry. 
He explored the basic structure responsible for the reality of the spectrum of 
non-Hermitian Hamiltonians, and established that all the PT-symmetric non- 
Hermitian Hamiltonians are pseudo-Hermitian. He also has shown that any 
diagonalizable operators with discrete spectra, is pseudo-Hermitian if and only 
if its eigenvalues are either real or grouped in complex-conjugate pairs (with 
the same multiplicity). Moreover this result has been generalized to the class 
all PT-symmetric standard Hamiltonians having M as their configuration space, 
and to the class of possibly nondiagonalizable Hamiltonians that admit a block- 
diagonalization with finite-dimensional diagonal blocks. In fact many of the 
later developments in the field are anticipated in the paper by Scholtz, Geyer 
and Hahne [18] (see also [19]). 

By definition [8] , a Hamiltonian if is called pseudo-Hermitian if it satisfies 
the relation 

ii f = rjHif 1 , (3) 

where 77 is a linear, Hermitian, and invertible operator. One can also express 
the definition ([3]) in the form if# = H, where H# = r}~ l H^f] is the 77-pseudo 
adjoint of if [5]. 

An interesting area where the pseudo-Hermiticity is applied is in the study 
of non-Hermitian two-level systems 114] . These simple Hamiltonian systems 
model accurately many physical systems in condensed mater, atomic physics, 
and quantum optics [20] - [26j . The latter field provides a beautiful imple- 
mentation of the coherent states formalism [28] - [31] . Rabi oscillations in the 
non-Hermitian system of a two-level atom in electromagnetic field have been 
recently examined in ref. [14j . 

In the preceding paper [32j . we have shown how the exact evolution and 
nonadiabatic Hannay's angle of Grassmannian classical mechanics of spin one 
half in a varying external magnetic field, is associates with the evolution of 
Grassmannian invariant-angle coherent states. 

In this paper we extend the fermionic coherent states approach [33] [34l [35l 
37J to two-level non-Hermitian Hamiltonians which are pseudo-Hermitian (p- 
Hermitian) . The underlying number system is Grassmann algebra }38 [ I39 | . The 
set of coherent states (CS) for pseudo-fermionic (shortly p-fermionic) system 
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turned out to consist of two subsets of states, which are bi-normalized and bi- 
over complete (shortly bi- normal CS). 

The paper is organized as follows. In section 2, we study a non-Hermitian 
two-level systems (a two-level atom interacting with electromagnetic field)and 
its pseudo-Hermitian properties. Then we introduce the creation and annihi- 
lation operators for the two-level p-Hermitian system with real energy spec- 
trum, such that it's Hamiltonian ascribes a form similar to that of the free 
harmonic oscillator: H = £l(b#b — 1/2), where b and b# are pseudo-fermionic 
(p-fermionic) lowering and raising operators. In section 3, we construct the 
p-fermionic CS as eigenstates of two annihilation operators, the eigenvalues be- 
ing complex Grassmann variables. The set of such eigenstates form a bi-normal 
and bi-overcomplete system. Then, in section 4, we study the time evolution 
of the constructed p-fermionic CS for the corresponding two-level p-Hermitian 
system. This evolution is shown to be temporally stable. The paper ends with 
concluding remarks. 



2 Two-level systems and pseudo-Hermitian 
fermions 

We consider a two-level atom interacting with an electromagnetic field. The 
general state of the two-level atomic system is 

m = c a (jt)\+)+c b (t)\-), 

where dp are the amplitudes of being in |=p). They are time dependent due 
to atom-field interaction. In the interaction picture (dipole interaction and 
phenomenologically described decay), and in the rotating wave approximation, 
the evolution of the system is described by the equation [24J [25] 



£ ( CM "<7. \ ( C a (t) \ 

l dt \ C' b (t) J- 2 \ * ~H b J \ C b (t) J ■ 

The real constants 7a , 7b are the decay rates for the upper and lower level 
respectively. The quantity uj characterizes the radiation-atom interaction matrix 
element between the levels ( lo* is the complex conjugate ). The basic vectors 
of the upper (lower ) level are | +) and | — ) . 

We remove the average effect of the decay terms by means of a nonunitary 
transformation in the state space, 

|V) - U{t)\i/>), U{t) = e rt , r - i( 7o + 7b ). (5) 

The probability amplitudes in the new representation are Ci(t) = exp(rt)C^(£), 
i = a,b, and the non-Hermitian Hamiltonian takes the following matrix form 
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where <J = ( 7o - 7fc )/2. 

The trace of H, eq.©, 



is vanishing, and the determinant of H is real, 



detH = (S 2 



! )/4. Therefore it is 77-pseudo-Hermitian (?7-p-Hermitian) [10] . 



Its matrix is a particular case of a more general 2x2 traceless matrix studied 
in [TT], where the complete biorthonormal system {\ipi), \<fri)} for H and the 
operator rj are explicitly constructed. We reproduce this system and r\ (up to 
certain common factors) in our specific notations. 

The eigenvalues Ei of H, i = 1,2, and the related complete biorthonormal 
system are given by (we consider the nondegenerate case of Ei 7^ 0) 



Ei 



2 



n 

~2 



(7) 



lV>i> 



'2Vt 



— lu* \/Q-\-iS 

M 



1*2 



1 I M 



'2Q 



(8) 



-LP* y/Q* —iS 

H 



id 



, (9) 



iS 



where f2 



U3\ 



For both real and complex eigenvalues (i.e. real and complex Q) the Hamil- 
tonian ^ satisfies the p-Hermiticity relation ^ with rj given by 



n = 



iSuj 



1 



(10) 



As noted in [T3] , the real eigenvalues correspond to the case where the dipolc 
interaction is large relative to the damping effect. In this case the ordinary Rabi 
frequency is replaced by the 'pseudo Rabi frequency' which in our notations are 
I a; |/2 and 0/2 correspondingly. 

In the case of = (i.e. \u}\ 2 = S 2 ), the amplitudes C' a ,C' b are given [M] by 
(1 — St) exp(— Ft) and (iuit/2) exp(— Tt) correspondingly, where T = ( 7a + 7fc)/4. 
Due to the exponential decay factor exp(— Tt), any divergence [27] does not 
occur in our system. 

In the case of \uj\ 2 < S 2 the eigenvalues of H are pure imaginary, but the 
Hamiltonian is still pseudo-Hermitian [14j . The PT-symmetry of our Hamilto- 
nian (6) is considered in p3] following the method of Bender, Brody and Jones 
[3] . The Parity operator P of the two- level system can be defined as [5] 



P 



1 

1 
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so that P 2 = 1, P = P- 1 . 

The generalized parity operator for the two-level systems is defined [T2J. US] 
as P = \(j)i){4>i \ — \(j> 2 ) (02 1 j which m °ur case results to 



P = 



"R 



v M 



Different definitions have been introduced by Mostafazadeh [TJ] and Ahmed 
[T3] for the antilinear time-reversal operator T. As in the paper [14], we use the 
representation introduced by Bender et al [3], namely, 

T = K 0l 

where Kq is the complex conjugation operator. One has Kq = 1, (PT) 2 = 1, 
and one finds that PT commutes with H: PKqHKq 1 P^ 1 = H. If one uses the 
generalized parity operator P = \<f>i){<j>i\ — \4 > 2)( < t } 2\> then one can take T = UKq, 
where U is a 2 x 2 unitary diagonal matrix, U22 = Uu = uj/\lj\. 

The generalized charge conjugation operator C for two-level system is given 
by the expression [TJ] C = \t(>i)(4>i \ — I "02 ) {^2 1 7 which in our case reads 

2n \-w -id J' n 

From the last equation we deduce that C commutes with the Hamiltonian 
H: [C, H] = 0. This invariance property eliminates negative inner products 

m- ' 

Furthermore, unless otherwise stated, we consider the small damping effect 
case of our system, i.e. lo 2 > S 2 , that is real (and if real it is positive). One 
can verify that, in this regime (i.e. \lo\ 2 > S 2 ) this 77, eq. (10), can be represented 
in terms of 1^) as 

*7 = o(l0iWil + l&Kfcl) = /W' (11) 

The operator r\ is positive definite. The construction i] + = l^iX^il + I^K^I 
was introduced by Mostafazadeh [IP] . 

Now, let us introduce the annihilation operator b associated to the Hamil- 
tonian H given in eq. ([6]), 




M 
M 



Its adjoint operator reads (fl is real) 




(12) 



(13) 



M 

and its 77-p-Hermitian adjoint b# , defined by 
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b* = r 1 - 1 b + r 1 , (14) 

takes the form 

b *-±( " M ^) (i5) 

Next we examine the properties of the operators b# and b. First of all is 
that b# and b realize a pseudo-Hermitian generalization of the fermion algebra, 
namely 

b 2 = b* 2 = 0, {b,b*} =bb* + b*b=l . (16) 

b& and b could be called the creation and annihilation operators of p-Hermitian 
fermions |40j . One can verify that they raise and lower the eigenvalues of H by 
a quantity £1 = 2E, i.e. they act on the states |^> 4 ) as follows 

6|^)=0, & ]-0 2 > = |-0x> > (17) 

b*\^ 2 )=0, b*\il) x ) = |V> 2 > , (18) 

The operator b annihilates the lowest eigenstates \tpi), and b# brings this state 
onto the upper eigenstates |-0 2 ) . 

Introducing the quadratic operator N = b#b, the p-fermionic number oper- 
ator, we find the following natural anticommutation relations 

{N,b} = b {N,b*}=b* . (19) 

In terms of the operators b and b# the Hamiltonian H is factorized (up to an 
additive C-number term) to a form, similar to that of the free (boson) harmonic 
oscillator, 

H = (b*b - ij . (20) 

Taking the Hermitian conjugate of both sides of (|2H)) we reaffirm the p- 
Hermiticity of H (according to definition ([3])): 

H + = nib+ribif 1 - -) 

— ^7777 1 (b + rjbrj^ 1 — — )r)rj _1 

The above relations confirm that b and b# are lowering and raising op- 
erators for the two-level p-Hermitian system (with real eigenvalues) and can be 
regarded as p-fermionic annihilation and creation operators. This is consistent 
with the limit 5 = 0, corresponding to a Hermitian Hamiltonian, when r\ = 1 



G 



and = b + , i.e. the p-Hermitian generalization of the fermion algebra reduces 
to the usual fermion algebra. Quantum system with Hamiltonian of the form 
(|2"0|) should be referred to as p-fermionic oscillator. 



3 Pseudo-fermionic coherent states 

Having introduced the p-fermion lowering and raising operators we now em- 
bark on the construction of the p-fermionic coherent states (CS) for our sys- 
tem described by p-Hermitian Hamiltonian H given in eq. ©■ We shall 
follow as closely as possible the scheme of fermionic CS developed in papers 
[331 [331 [351 [37] , generalizing it to the p-fermion case. 

For the reader convenience we begin with a brief reminder of the properties of 
complex Grassmann variables [331 EH [36l [37] , denoted here as £ and £* . 

The complex Grassmannian variables ^ and their complex conjugates £* 
satisfy the anticommutation relations: 

tei+Z& = °> (21) 
{£,^}=0, {£*,£*} =0. (22) 
The £j 's anticommute with b and b# , 

{Zi,b} = 0, {£*,&} = 0, {&,&#} = (), (23) 
and have the following properties 

£IV>i> = h/>i>£, ei^ 2 > = -|-0 2 ) , (24) 
^i) = \4>i)L £l*a> = -|&>£ ■ (25) 

The pseudo-Hermitian conjugation reverses the order of all fermionic quantities, 
both the operators and the Grassmann variables: 

(b*Z i + eib) # =e i b + b%. (26) 

The Grassmann integration and differentiation over the complex Grassmann 
variables are given by 

#1 = 0, J dtf = 1, J dC 1 = 0, j dCC = 1 , (27) 
d d _ d d j . 

^ 1 = < ^ =1 > w 1 = ' w c = 1 - (28) 

The Grassmann integration of any function is equivalent to the left differentia- 
tion 

/de/(0 = ^/(£). (29) 
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We define the displacement operators D(£) for any set of complex Grass- 
mannian variables £ in the following way: 



D(t)=exp(b*£-?b) 

= l + b*£-Cb+ [b*b 




(30) 



(31) 



The pseudo-Hermitian adjoint D# is given by 



D#(0 =cxp (tb-b*£) 

= l+Cb-b*£ + (b#b-V\z*£. 



(32) 



(33) 




Using the explicate formulas of D and D# , and the anticommutation rela- 
tions between operators b, b# and Grassmann variable £ we establish that D(£) 
are pseudo-unitary: D*(£)D(£) = 1 = D(£)D*(£). 

We now define the p-fermionic CS |£) as eigenstates of the annihilation op- 
erator b, 



The eigenvalue £ is a complex Grassmannian variable. 

The Hermitian adjoint of the CS (the bra- vector) is (£| and it is left eigenstate 

of b*, (£\b' = In order to meet the alternative relation n{£,\b^ = t;(£|£* 

one has to define = (|£)) # := (dis- 
similarly to the cases of Glauber bosonic CS [28] and of fermionic CS [33] our 

p-fermion eigenstates |£) can be constructed from the lowest (ground) eigenstate 

\ipi) of the Hamiltonian H, acting on it by the pseudo-unitary operator -D(£): 



By using the formula pip for the displacement operator, we may write the state 
|£) in the form 



MO = £10- 



(34) 



10 = D(0 |^) . 



(35) 




(36) 



The Hermitian adjoint of the CS is 



<£l = Wi 1 1^(0 



(37) 




(38) 
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and the inner product is 

(tit) = (i>M + «V> 2 |^ 2 ) - (V*»£*£ - 2fljn(e<^ 1 |V 2 » . (39) 

so that the |£) are not normalized. 

Now we have to examine for (over)completeness the set One can 

straightforwardly check (using the rules (f2"Tj) - (HZ])) that neither the inte- 
gral (against the measure d£,*d^) of the Hermitian nor the integral of 
p-Hcrmitian |£)tj(£| (unnormalized) projectors result in the identity operator: 

J «|£><£| ¥= 1, J v(t\ + 1 • (40) 

The way out of this impasse is suggesting by the known transition from 
'orthonormal system' of eigenstates of Hermitian H to the 'biorthonormal sys- 
tem' of states of p-Hermitian H. With this idea in mind we introduce another 
continuous family of states namely the eigenstates |£) of the operator b, that 
annihilates the dual state \(j>i), 

b\0=m, h = iPrT x . (41) 

Operator b is nilpotent, b — and anticommutes with w . Representing w = 
rjWri^ 1 we see that b^ is ry'-p-Hermitian adjoint to b, rj = rj^ 1 . Denoting this 
pseudo-conjugation by we obtain the pair of p-fermionic operators b and 

bb*' + b*'b = l, b 2 = {b*') 2 =0. (42) 

In view of the p-fermionic algebra (|42|l we introduce new displacement op- 
erators 

D(0 = exp(&#'£ - C~b), D*'(£)bD(0 = b + £, 

and construct eigenstates of b according to the above described scheme (see eqs. 
(ESI), ©), 

\0 = D(0\4>i) (43) 
= exp (I0i> - €l^a» - (44) 

The scalar product between takes the form 

m = (4>M + «^2> - - 2ilm(£<<^ 2 », (45) 

while 

or more generally, 

<<eio = <^|Dt(o5(cM) = rc + j(2 - ro(2 - c*c). (46) 
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By means of the two type of states |£) and |£) the resolution of the identity 
is realized in the following way, 

i = J drdeiexei = / «li><e (47) 

The equations (|47|) can be easily verified using the expansions of |£) and |£) in 
terms of \4>i) and (eqs. (136]) and ((44])) and the rules of permutation and 
integration - (j2T)) . 

We have obtained that the system of one-mode p-fermionic CS consists of 
two subsets, namely and l n view of ([45]) and (|46[) this continuous 

system should be called bi-normalized and bi-overcomplete, or shortly system 
of bi-normal CS. Similarly the two sets of pseudo- unitary operators D(£), -D(£) 
should be called bi-unitary: 

£(0^(0 = 1 = ^(0^(0- 

Note that £>(0 is 77-pseudo-unitary, while D(0 is f/-pseudo-unitary with rj' = 

4 Time evolution of p-fermionic coherent states 

A given parametric set of states is said to be realizable for a physical system 
if the time evolution \ip; t) of any initial state from the set, governed by 
the Hamiltonian, leaves the state in the set [42] • In other word \ip',t), for any 
t, obeys the defining criteria of the set. In such a case one shortly says that 
the time evolution (of the parametric set of states) is stable [32] • In Hermitian 
mechanics this means that the time dependence of the states is included, up to 
a phase factor, in the state parameters. For example the time evolution \a; t) of 
Glauber CS |a) [28] is stable with respect to the harmonic oscillator evolution 
operator exp(-iHt), H = uj(a^a + 1/2): 

\a;t) = e' iuJt/2 \a(t)), a\a;t) = a(t)\a(t)), a(t) = ae~ iu}t . (48) 

In the case of our p-fermionic CS the set' parameter is the complex 

Grassmann variable £, the eigenvalue of the p-fermionic lowering operators b or 
b. The time evolution is stable if the evolved states |£;t) and |£;t) remain 
eigenstates of the operators b and b respectively, 

This implies that the time evolved CS {|£;t) and \£;t)} should form bi-normal 
and bi-overcomlete system. 

Let us first consider the time evolution of an initial CS Clearly we have 
|0 i) = &q?(-iHt)\£), |0 0) = |0- Using the form ((36j) of |0 and the facts that 
\tjj 1 2 ) are eigenstates of H (with eigenvalues £1,2) we get 
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16 t) = e-^ (l - |lM - e-^^|^ 2 > • 

Taking into account that £1 = —£7 and E2 = E we put £(t) 
rewrite the last equation in the form 

= e iEt \m- (51) 

which manifests the stability of the time evolution of CS Note that the 
overall time dependent factor exp(iEt) is a phase factor since Ei are real. 

In a similar manner we establish, that the time evolution |£; t) of an initial 

is stable (remains eigenstate of b): 



(50) 
£ and 



IC;*> 



„iEt 



£(*)*£(<) ) i^i> - ec*)i^ 2 



exp(-vEit) ( 1 - ) I'M " exp(-iE 2 i)£|0 2 



|€;t)=exp(-ifl+t)|0 

~ 2 

= exp(zM) ^1 - ±£(tyt{t)\ 1^) - £(W 2 >) = exp(iEtMt)) • (52) 

The results (|5"Tj) and (j52"|) reveal the bi-normality and bi-overcompleteness of the 
family of time evolved states {|£;t), \£',t)} of the p-fermionic oscillator system 
P0|) : one has = 1, and 

1 = J dCd^ ; t)(M\ = J de*deief*><*; • (53) 

We observe that here the time evolved states |£;t) and |£;t) differ from CS 
|£(t)) and m phase factors only. In more general cases the overall factors 

Af(t) and jV'(t) ascribed in the stable evolution of bi- normal and bi-overcomplete 
system of states could not be phase factors, but their product should equal unity, 

M*{t)M{t) = l. 

Finally we have to note that a complementary bi-normal and bi-overcomplete 
system of states can be constructed, in a symmetrical manner using the opera- 
tors 6* and b#, that annihilate the 'upper-level states' 1-02) and |<fe)- 



5 Concluding Remarks 

In this paper, we have generalized the fermionic coherent states (CS) for 
two-level systems described by pseudo-Hermitian Hamiltonian with real spec- 
trum. Unlike the standard bosonic and fermionic cases the system of pseudo- 
fermionic (p-fermionic) CS consists of two subsets, which are bi-normalized and 
bi-overcomplete. In this sense the system of p-fermionic CS can be regarded as 
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a continuous analogue of the bi-orthonormal system of discrete eigenstates of 
p-Hermitian H. The two subsets are built up as eigenstates of the p-fermion 
annihilation operators b and b — r/h]^ 1 , where r\ is the Hermitian operator that 
ensures the p-Hermiticity of the Hamiltonian, H — r^ x }i^r\. In terms of b and 
5# = ^-ifrt^ the Hamiltonian is factorized to the form of p-fermionic oscillator, 
eq. (HOD. 

The evolution of the p-fermionic CS governed by the p-Hermitian two-level 
Hamiltonian © is shown to be time stable - the evolved states remain eigen- 
states of the p-fermionic annihilation operators, preserving the bi- normality and 
bi-overcompleteness of the system at later time. In the Hermitian limit of ?y = 1 
(that is 8 = in (JBJ) our p-fermionic CS and all related formulas recover stan- 
dard fermionic CS of refs. [33l [34] . Time evolution of fermionic CS for a Pauli 
spin in a slowly varying magnetic field was examined by Abe [35] (see also the 
comment [36]). 
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